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SDMKARY 


On the basis of a recently developed theory for sveptback wings 
at supersonic velocities, equations are derived for the vave drag of 
sveptback tapered wings with thin symmetrical double-wedge sections 
at zero lift. Calculations of section nave -drag distributions and 
wing wave <jtrag are presented for families of tapered plan forms. 

Distributions of section wave drag along the span of tapered 
wings are, in general, very similar in shape to those of untapered 
plan forms. ‘Far e given taper ratio and aspect ratio, an appreciable 
reduction in' wing wave-drag coefficient with increased sveepback is 
noted for the entire range of Mach number considered. For a given 
sweep and taper ratio, higher aspect ratios reduce the ving wave- 
drag coefficient at substantially subcritical supersonic Mach numbers. 
At Mach numbers approaching the critical value, that is, a value equal 
to the secant of the' sweepback angle, the plan forms of low aspect - 
ratio have lower drag coefficients. 

Calculations far wings of equal root bending stress (and hence 
different aspect ratio) indicate that tapering the wing reduces the 
wing wave -drag coefficient at Mach numbers considerably less than 
the critical value but increases the drag coefficient at Mach numbers 
near the critical values. Comparisons on the basis of constant aspect 
ratio, however, indicate an increase of the wing vave -drag coefficient 
with taper at Mach numbers considerably less than the critical value 
and a decrease of the drag coefficient vith taper at Mach numbers 
near the critical value. 


INTRODUCTION 

Recent developments in airfoil theory for s\tpersonic speeds 
(references 1 and £) indicate pronounced favorable effects of sweep- 
back on the wave drag. In reference 1, a method is developed for 
calculating pressure drag at supersonic speeds for sueptback airfoils 
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having thin sections at zero lift. Reference 3 applies this method 
to calculate the supersonic wave drag for a series of untapered wings 
with symmetrical "biconvex airfoil sections. 

The present paper applies the method of reference 1 to derive 
the generalized equations for the section wave drag and. wing wave 
drag of sweptback tapered wings with thin symmetrical double -wedge 
sections at zero lift. Section wave-drag distributions and wing 
wave-drag calculations are presented for specific families of tapered 
plan forms. The airfoil sections and wing tips are chosen parallel 
to the direction of flight. The angle of sveepback is referred to 
that of the line of maximum thickness, and the range of Mach number 
considered is between 1 and the critical value corresponding to the 
condition vhere the Mach lines are parallel to the maximum-thickness 
line, that is, to a Mach number equal to the secant of the sweepback 
angle . 


SYMBOLS 


x, y, z 

V 

P 

Ap 

<1 

<P 

M 


cartesian coordinates 
velocity in flight direction 
density of air 
pressure increment 


dynamic pressure 



disturbance -velocity potential 
Mach number 


p = \|m 2 1 

dz/ax slope of airfoil surface 

a root semichord, measured in flight direction 

c chord length at spanwise station y, measured in flight 

direction 

'S. 

maximum thickness of section at Bpanwise station y 


t 
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A 

angle of 

®0 

slope of 

m l 

slope of 

“2 

slope of 


sweep of the line of maximum thickness, degrees 
line of maximum thickness (cot A) 


wing leading edge 


wing trailing edge 


2*1 ~ “ 0 ) 


d= 

“0 


c tip 




Oft 


tip 


CD 


epan of wing 


■wing area 

aspect ratio 

taper ratio, ratio of tip chord to root chord 

section vave-drag coefficient at spanviBe station y 
exclusive of tip effect 

increment in section -wave-drag coefficient at spanvise 
station y due to tip 

section wave-drag coefficient at span-wise station y 

( c w + Ca t iP ) 

wing vave-drag coefficient exclusive of tip effect 
increment in wing wave-drag 
wing wave- drag coefficient 


coefficient due to tip 

( C D„ + C B t ip) 



Subscript s refers to conditions at root 
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ANAIXSIS 


The analysis is based on supersonic thin-airfoil theory and on 
the assumptions of small disturbances and a constant velocity of 
sound throughout the fluid. These assumptions lead to the linearized 
equation for the velocity potential cp (reference 4) 



where M is the Mach number of the flow and the derivatives are 
taken with respect to the variables x, y, and z of the cartesian- 
coordinate system. It should be noted that the linearized theory 
is not expected to be applicable near Mach number unity. On the 
basis of this linearized theory, a solution for a uniform swept- 
back line of sources in the pressure field is derived in reference 1. 

The pressure field associated with this solution corresponds to 
that over an airfoil of wedge section. The pressure coefficient 
Ap/q at a spanwise station y and point x along the wedge is 

££ = £ 45. — cosh- 1 x , - (2) 

1 {1 - eV e|y-»ix| 

where m^ is the slope of the leading edge of the wing, dz/cbe is 
the tangent of the half-wedge angle (approx, equal to half-wedge 

angle since the angle is small), (3 = yM 2 - 1 and the origin of 
the line source is taken at (0,0). 

The distribution of pressure over sweptback wings of desired 
plan form and profile is obtained by superposition of wedge— type 
solutions. In order to satisfy the boundary conditions over the 
surface of a tapered wing of symmetrical double-wedge section, semi- 
infinite line sources are placed at the leading and trailing edge of 
the wing and a semi-infinite line sink of twice the strength is placed 
along the line of maximum thickness so that all three lines intersect 
at one point. At the tip where the wing is cut off in the flight 
direction, a reversed distribution of these lines of sinks and sources 
are p3.aced so as to cancel exactly all effects of the original distribu- 
tion farther spanwise than the tip. Figure 1 shows the distributions 
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of BlnlcR and sources for a tapered ving and identifies the system 
of axes and the symbols associated 'with the derivation of the drag 
equations. 

The disturbances caused by the elementary line sources and 
sinks are limited to the regions enclosed by their Mach cones. 

Figure 2 shows the Mach line configuration for the tapered-ving 
plan form and indicates the regions of the wing affected by each 
line source and sink. For purposes of simplification the tapered 
vings considered vere restricted to those -with no tip effects other 
than the effects each tip exerts on its am half of the -wing. For 
a ving of taper ratio 0, no tip effects need be considered since 
the Mach lines originating at the tip do not enclose any part of the 
wing. 


The pressure coefficients obtained from superimposing solutions 
of idle type shorn in equation (2) are converted into drag coeffi- 
cients by the following relations: 

For section drag at a spanwise station y 


c*c 


p Trailing edge 

Ap dz 
q dx 

Leading edge 


dx 


( 3 ) 


where 


y/^m-L - mg) + 

C 5* . 

m^ 


is the chord length at y, and the integration is performed along 
the chord parallel to the flight direction. 
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The ving wave-drag coefficient is obtained by integrating the 
section drag along the span and dividing the result by the ving 
area . 


C 


D 


2 


n Tip 


s I 

UEoot 


c d c 



ITip HT.E. 
iJRoot UL.E. 


Ap dz . . 

— — dx dy 
q cbc 


w 


where S is the ving area, and the integration with respect to y 
is performed along the span. 


DERIVATION OF GEEERALIZED EQUATIONS 


By superposition of wedge-type solutions (equation (2)), the 
pressure field is obtained for a tapered wing with leading edge, 
trailing edge, and line of maximum thickness sweptback. The drag 
equations are derived for half of the wing 3ince the drag is distri- 
buted symmetrically over both halves. The induced effects of the 
opposite half -ving are represented by the conjugate terms in the 
integrands of the drag integrals. 

For a symmetrical double -wedge profile, 


dzl 

dx 


t 

c 


where t/c is the section thickness ratio. The generalized equa- 
tion, exclusive of tip effects, for the wing wave drag is obtained 
as follows: (See fig. 3 for information pertinent to integration 

limits. ) 
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TtSC- 


% 


P^O 


8(t/c) S 4(t/c) 2 J 0 


Cfl C fly 


m l 






i UO J 


y-am 1 


A dx dy 


m. 


iamo 


JO J 


y-s-asig 

~ 


A dx dy ! 


y_ 

mo 


2no 




«% 

1-prin 


UO JPy 


®0 


pamo 


B dx dy - 


00 


y+amg 

mg 


B dx dy 


mo 


ndaiQ 


auh 


mo 


B dx dy 


y-a% 


l-pm^ 


8 3 33q jT-r OIEg 


“2 


\/i - p 2 v 


l-p^O 


mg 


C dx dy 


UO JPy+a 


nSSL nJL 


W*! 

SWq 


•p^o 


“0 


C dx dy + 


uPy-'-s 


Idmo 

2am.^ 


n — 

mo 

[y-^ 

m-i 


C dx dy 


^dmo 


amo 


y+amo 


mg 


C dx dy 


^I-Pjho l ffio 


(5) 
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vhere A, B, and C refer to the pressures resulting from the leading 
line sources, line sinks, and trailing line sources, respectively. 


2 o 

x+a+m-ipy , x + a - m, p y 

A = cosh x - + cosh x i 


P |y + ffi] 


x + am. 


Pjy - m-jX - am .2 


_ x x + mop 2 y , x - moP 2 y 

B = cosh 4 cosh 


p|y + ihqX j 


P|y ” mox| 


-1 x - a + mpP 2 y -i x - a - MoP 2 j 

C = cosh x — — + cosh 1 1— : 


P jy + BigX - anig j 


p y - n^x + amg 


The limiting case (taper ratio 0) is obtained by letting d = — 

“0 ‘ m l 

and the ving of constant chord (taper ratio 1.0) is obtained by 
equating m, = icq = ru. The integrations in equation (5) are per- 
formed and the resulting formulas for the section vave drag and the 
ving vave drag for the complete range of conventional taper (0 S taper 
ratio 5 1.0) are presented in appendix A. 

It vas stated previously that the tapered vings considered have 
no tip effects other than those each tip exerts on its ovn half of 
the ving. This implies that the Mach lines from one tip do not 
enclose any part of the opposite half -ving. This condition is 
expressed mathematically as follows: 

For Puig < 1 


Aspect ratio = 


2dtriQ 


2m.mQ 


a(l + X) (l + l)(Pm 1 mQ + idq - m-^) 


and for Pm^ > 1 


2dm0 km. 

Aspect ratio = > 

a(l + X) " (1 + M(1 + Pm-J 


(6a) 


(6b) 
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■where X is the taper ratio f ■ ^ or ^’ \ it can he seen from 

\E 00 t chord/ 

equations (6) that this simplification does not materially limit the 
range of Mach number that may he considered. Far small taper ratios 
this limiting effect is negligible and for taper ratio 0 there is no 
limitation whatsoever since equations (6) reduce to expressions that 
are alwayB valid. 

The wave-drag contribution of the tip is (see fig. 3) 


^ip « c cdy 

8(t/c) 2 Mt/c) 2 am 0 ia 2( 1 ''' |3lll l)~ gam l ln g dtlp 

m 1 (l+pm 2 ) 


\jl - p 2 ^ 2 J aa 0 g ( 1+ P a l)" am l 2D 0 dm 0 (l+Pm 1 )"m 1 (a-Hpy) 
|_ m l(l + P B, o) 1 m l 


dEQ 2 ( l+Pm^ )-amT_mo ? ' ^ 

m^i+pmo) 1 2q 


dm3 2 (l-!-Pmi )-amimQ y+amg 


m^l+pmQj 


I dmcp2(l+Pmi)-2amimg | toiQ ( l+pm^ ) -m^ ( a +py ) 
m^( l+pBig ) m l 


D dx dy 


dry, 1 

a°o ” “2 


E dx dy 
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where D and E refer to the pressures resulting from the leading 
line sink and line source, respectively . 


D = cosh 


m 1 (x + a) - d m Q - m 1 2 p 2 (y - dmp) 
Pm-Jy - m x (x + a)| 


E = cosh 


x ~ d " p 2 m o( y ~ dm o) 

P|y - no* | 


The Mach cone from the trailing line sink at the tip does not enclose 
any part of the wing and, hence, has no effect on the we drag. 

Equation ( 7 ) is solved for section wave drag and wing wave drag 
for the complete range of taper and the results are presented in 
appendix B. The total wave-drag coefficients are then obtained by 
the following relations: 


c a 3 c a* + c ati P 

c D » c Doo + c Dup 


(8) 


It is found that Ct) , is identically eQ.ua 1 to zero for all 

t xp 

cases satisfying the aspect -ratio limitations expressed in equa- 
tions ( 6 ) and, hence, C D = C D for the tapered wings considered. 


The conditions imposed in equations ( 6 ), although not materially 
limiting the range of Mach number for tapered wings, do limit to a 
certain extent the range of Mach number for low-aspect-ratio wings 
of constant chord. Equation ( 6 a) for this case reduces to 


Aspect ratio 



since m^ = iBq = 
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For untapered v3ngs of aspect ratio 2, 1, and 0-5, the lowest 
Mach numbers that can be considered without talcing into account 
additional tip effects ere 1 . 118 , 1.414, and 2 . 236 , respectively* 

It is desirable, therefore, to take into consideration for untapered 
plan forms the induced effects of the opposite tip when the Mach 
lines from one tip enclose parts of the opposite half-wing. Figure 4 
shows the Mach line configurations for these induced effects, and the 
drag equations are derived in appendix C. The ■'.ring wave-drag coef- 
ficient is then obtained from equation (8) -where Oj^_^ for these 

cases includes the effects Induced by the opposite tip. 


ESSUUS AND DISCUSSION 


Calculations were made for families of tapered plan forms, each 
family characterized by a constant sveepback of the maximum-thickness 
line. The plan forms were obtained by considering the moment of the 
area about the root chord divided by the cube of the root chord to 
be constant for any given family. The aspect ratio varies with 
taper ratio because of this area -moment parameter. 


For a constant thickness rat?o the parameter, area moment 
divided by the product of the root chord and the square of the root 
thickness, is also constant. This condition is intended to imply 
that to a first approximation the root bending stress is the same 
for all members of any family having the 3ame thickness ratio. A 

«</i ^ r» ^v i, 4 r>4.4 a v«a *1 -,r jwP A -PnyTWfl I-1 ■f '. p n yflyf Q — 

iu V ^ j aJaxx-i^ ui iJ"P"rou j7.£2iri x ttttbxs a lit- c tyjptyw u xa xrXv vcjx x-a 

tion with taper ratio Is shewn in figure 5* 


Section wave drag .- Section wave -dreg distributions for wings 
of taper ratio Q 0.5, and 1.0 are presented in figures 6 to 10 for a 
Mach number of 1.414 and streepback of 60°. The distributions of 
section wave drag of tapered wings are, in general, very similar to 
those of untapered plan forms. As a point of interest, the induced 
effecte of the opposite half— wing and the tip-effect distribution 
are shown in figure 10 as separate curves. The total section wave- 
drag distribution is then obtained by adding the tip drag curve to 
the solid-line curve. The tip effect is placed correctly as shown 
for a wing of aspect ratio 1.0; for a wing of aspect ratio 2, this 
tip drag distribution should be shifted 1 semichord to the right. 

It is seen by reference to figure 5 that figures 6, 8, and ID 


/ -P* r. 
\“b 


1 \ ^ a+v'-t’Kll+’i Ar»o -PrVr* ^OW! ItT 

_ f LAX <3 oouwi Uli nu i v uj. u £5 u-J. u WJ. XMM.WX A V* -*» ».*»*»*.> -V 


of wings and that figures 7, 9, and 10 (fig. 10, A = 2) are for 
another family of wrings whose aspect ratios are twice as large, 
respectively. 



12 


NACA EM No. L7E23a 


It is interesting to note at this point that for a given Mach 
number the section wave-drag coefficient at the root is a function 
of the sweep of the maximum-thickness line only; the terms involving 
leading-edge sweep adding up to zero. (See section drag equation 
in appendix A for y = 0.) 

Wing wave dr ag.- Typical variations of wing ■vave-drag coeffi- 
cient with Mach number for wings of taper ratio 0 and taper ratio 1.0 
of the same family are shown in figures 11, 12, and 13 for 50°, 60°, 
and 70° sweepback, respectively. At some Mach n umb er between 1.0 
and the critical value (M cr itical = s©c A), the drag curve for the 
tapered wing has a discontinuous slope. This discontinuity occurs 
at that Mach number corresponding to the condition where the rear 
Mach line crosses the trailing edge of the wing, that is, where 


mg mjmQ 


In this region and near the critical Mach number 



the 


theory is not expected to be applicable because the assumption of 
small disturbances is violated, but the results are presented in 
order to give a more complete picture of the linearized theory. 


It is seen from figures 11 to 13 that taper reduces the wing 
wave-drag coefficient at Mach numbers substantially below the 
critical value but increases the drag coefficient at Mach numbers 
approaching the critical value. This trend is similar to the one 
shown by the effect of high aspect ratio on the wave-drag coeffi- 
cient of wings for a given taper ratio. It must be remembered 
that for the families of tapered wings considered in these calcula- 
tions, however, the wings with greater taper have higher aspect 
ratios and, hence, the effects of aspect ratio as well as taper 
are included in this trend. 


Variations of wing wave-drag coefficient with taper ratio for 
different sweepback angles at a Mach number of 1.2 are shown in 
figure 14. The untapered wing for this family has an aspect ratio 
of 1.0 and the variations of aspect ratio with taper ratio are • 
presented in tabular form in the figure. For a given sweep angle, 
the wing of taper ratio 0 has the lowest drag coefficient and the 
untapered wing the highest. As the Mach number approaches the 

critical value ( (3 = — ), this trend would reverse itself and the 

\ V 
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untapered wing would have the lowest drag coefficient as can "be 
seen by reference to figures 11 to 13 . It is also evident frcm 
figure 14 that for a given taper ratio and aspect ratio, an 
appreciable reduction in wing wave-drag coefficient is accomplished 
with increased sweepback. 

Figure 15 presents variations of wing wave -drag coefficient 
with taper ratio for three families of wings based on untapered plan 
forms of aspect ratio 0-5, 1, and 2, respectively . The results are 
presented far 60° sweepback and a Mach number of 1.4lh. Pertinent 
details of the wings are presented in tabular farm in the figure to 
facilitate interpretation of the plotted curves. The aforementioned 
trend of reduction in wing wave -drag coefficient associated with 
high aspect ratios at Mach numbers substantially below the critical 
Mach number for a given taper ratio is clearly seen in this figure. 
By choosing points along these curves corresponding to wings of the 
same aspect ratio, it is seen that for a constant aspect ratio 
tapering the wing increases the wing wave -drag coefficient. By a 
similar procedure it can be shown that for wings of constant aspect 
ratic* taper reduces the wing wave-drag coefficient at Mach numbers 
near the critical value. The increase in aspect ratio with taper 
ratio defined by the area -moment parameter thus has the effect of 
offsetting the adverse effects of taper at Idle lower Mach numbers. 


CONCLUSIONS 


1. Distributions of section wave drag along the span of tapered 
wings are. In general, very similar in shape to those of untapered 
plan forms. 

2. The section wave-drag coefficient at the root is a function 
of the Mach number and the sweep of the maximum-thickness line and 
is independent of taper. 

3- The increment in wing wave-drag coefficient caused by the 
tip is identically equal to zero for all tapered and untapered wings 
for which the Mach lines from one tip do not enclose any part of 
the opposite half -wing. 

For wings of equal root bending stress, taper reduces the 
wing wave-drag coefficient at Mach numbers considerably lees than 
the critical value - that is, a value equal to the secant of the 
sweepback angle - but increases the drag coefficient at Mach numbers 
near the critical value. 
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5* For wings of constant aspect ratio, taper increases the wing 
wave-drag coefficient at Mach numbers considerably below the critical 
value and decreases the wing wave-drag coefficient at Mach numbers 
near the critical value. 

6. For a given taper ratio and aspect ratio, an appreciable 
reduction in wing wave -drag coefficient with increased sweepback is 
noted for the entire range of Mach number considered. 

7« For a given sweep and taper ratio, higher aspect ratiqs 
reduce the wing wave-drag coefficient at substantially subcritical 
Mach numbers. At Mach numbers approaching the critical value, the 
plan forms of low aspect ratio have lower drag coefficients. 

The generalized equations presented in the appendixes may be 
used to calculate the subcritical supersonic wave drag at zero lift 
for any conventionally tapered or untapered wing with symmetrical 
double-wedge airfoil sections and with leading edge, trailing edge, 
and line of maximum thickness sweptback. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va . , 
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APPERDU A 


EVALUATION OP EQUATION (5) FOR SECTION WAVE IRAQ AND WHO 


WAVE DRAG EXCLUSIVE OF TIP EFFECTS 


N] 


Section Drag for 0 < Taper Patio < 1 


0 < y < 


1 - pn^ 


= A + B 


ami an^ 

< 7 < — 

1 - Ps^ *“ 1 - P»fe 


= A + B + C 


afflo ^ 2am l 

— < y < 

1 - Pn$ - 1 - Pmj_ 


A + B + C + B 


2eaar an^cFi 

< y < 

1 - pm^ ~ hiq - mi 


r- 

2 J j^O 

1 - p^nu 2 L 


+ mi ) + aswh _i y(l + no^iP 2 ) + “b 
— — cosh " ' • " 

no P(y(nb + m i) + an i D bj 


y cosh 


_l 1+ pV 2 yfnh - fflo) | + aqiPb cosh -I j(l ~ ^i^ 2 ) 1 


2?^ 


P jy(®o * m i) " ^l^bl 
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y ( m 2 + m i) + SamjBu , y(l + m-imoP 2 ) + 2amp 
— cosh — —A — 

^2 P |y(»l + mg) + 2am 1 m 2 J 


y( m l " “ 2 ) + 2 am i m 2 _] y(l - mnmoP 2 ) + 2auu 
— — cosh ' — - ^ 


21-12 


Pjyfmg - 2 ^) - 2am 1 m 2 | 


i 


2 y , -1 1 + ®2^P 2 

— cosh 


1 - p 2 m 2 2 


2PB12 


I 1 “ P 2ffl o 2 


Jy(mo + mg) + anting _ x y(l + mQiiigP 2 ) + auu 

cosh ^ 




p|y(mo + mg) + anyogj 


- % coeh -1 1 - >a 0 gpg /( m 0' m s)^^ 008h -l y(l-^ g ) ta% , I 
2Pm 0 “2 P h r ( m o ■ “o) “ aD b m 2 1 1 
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B * 




) r -| -! + ” o ” i » 2 ) ' 

0 ^ + - "H COBh 


+ [ 7(120 - a ^) - aryaj 


cosh 


-1 


y(l - a^ajP 2 ) - an^ 

0|y( m l - * 0 ) + a^ifflo 



c = 


“o’ 


if ~ P% £ 


] r , v n -1 yC 1 + “ o ^ P 2 ) " 

<[H*o + “ 2 ) - am 2 ffl o C08h ~n ; 

| L • p|y(n^ + mg) - a 


°“0 
ajn 2®0 | 


+ [y(m^ - mo) - amgSJo"! coBh' 


. x y(i ~ “o^P 2 ) ~ Wp 
P|y(“o " “2) + SJB 2 E b | 


and 


= — -— Jr ™, - M ^ 1 «** jfeljatfha . 

*i l fi - p 2 ^ 2 [ L - Piy( m i ■ “2) + 


- jy(m-L + nig) ■ cosh 


■1 y(* + “l^P 2 ) “ 2am l 


x AJ Bly/n-. + nu \ - Prsm~m-\ 

f \ J- <Z. r 


1 

* L " a | ' 
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Wing Drag for 0 < Taper Ratio < 1 


8(t/c y 


= A 


for 0 ^ dniQ ^ 


1 - pm-L 


= A + B 


ami _ , < “b 


1 - pm- 


_<dmo< — 


pxno 


= A + B + C 


-2»_ < - < Jf3_ 

1 - Pidq ^ 1 - pm x 


= A + B + C-t-D 


2ami <a^< ; aVl 


1 - ^ “o • m l 


•where 




2 Jm 0 -' m l a “ , d(l - p 2 ^ 2 ) + a .-.l+rn^p 2 

A <— cosh 1 2d 2 m C) <:j cosh 1 


Jl-P 2 ^ 2 |mo 2 -% 2 


a 


PmQ 


2pmo 


m 0 2 [m 2 (d - a) - d^] 2 ^ dm^l - n^p 2 ) + aim, 

21112(102 - mo) PmQjmgCd - a) - dmQ j 

idq 2 JdroQ + iQgCd + a) j 2 dm^l + m^P 2 ) + sa 2 1 

+ ; ; — cosh" 1 > 


2m 2 ( m 0 + “2 ) 


PmojdmQ + ingCd -h a ) J 


jmojdmo + mi(a + d)"| 2 _ 1 d(l + P^Mij + a 


cosh 


\Jl - P 2 m x 2 (_ 2 (“o + m l) 


+ mi(d + a)J 
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d 2 !^ 2 l + m^ 2 P 2 mg fcte^ - i&i(d4- a)] 2 , d^l - + a 

cosh — --- cosh~ x 

2 2 Pm^ 2^-21^^ Pjdn^-m 1 (d+a)| 


+ mg ) + Sam-jingj 2 dm^l + p^n^rog ) + 2amg 

Ifflig^ + mg) ppm^s^ * mg) + 2am 1 a^] 

JjdmQ^mg - m^) - Sam-^mgJ 2 <3 iiiq(i - rnjiBgp 2 ) + 2amg I 

lnng(mg - m x ) - a^) - 2am 1 m 2 | j 


* .S' 2 .-1 1 + ^g 2 


-a 2 *, 2 ooah' 1 __2J 

^ - p V L 


Ua'Sng^m-j^ ^ too(l “ nig 2 ? 2 ) + 2amg 


mg 2 - m x 2 


2pam 2 e 


2a 2 mg , ^o(l " EJg 2 ? 2 ) + 


+ -■ - ■ - ■ - cosh 

V 2 " *b 2 


W 


cosh” 1 -i- + 


(m 0 +m 1 )^.-p £ 


Pm l (3m 1 -nfe)\Jl-P^ 2 


cosh 


II. T cosh” 1 

l 3 ®! " “oH^O + m l)£ ” >V PD, ° 


*When 1 - mg 2 p2 ie negative neglect term marked vith asterisk 
(-d®nu2 cosh” 1 ---t— ^., - P- ) in values for A and use the relation- 

' l x 

Ship cos X x = -i cosh -x x for all terms involving — ■ — 

Jl - mp2p 2 

-I rto+4 />«« a4> ^ 


multiplication factor. 
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B = 


a2m i^\l 1 " “o 2 ^ 2 -i (ta o( 1 " m i 2 P 2 ) " 

cosh 


2mp J 

- P 2 eiq 2 ^m^ ” ~ 

w 


aPm-, 2 


ffiofamo + m^(d - a)l 2 , «^o(l + “o m ]P 2 ) " am l 

-Alls = — cosh 


?m 1 (Di 0 + m 1 ) 


Pjtiq Id^niQ + m 1 ) - am-L 


C = 


rnQ [!!]_( d + a) - dmol 2 dno(l - moiJiQp 2 ) - an^ 

— — cosh — — ■ 

2mi(mi - Mq ) 

2 xdq 1 ° 2 t„_ 2 - 


Pm 0 |d(m 1 - joq) ■■■ am 1 


a y mo 2 mp\jl - mg 2 P 2 d(l - mo 2 P 2 ) - a 

' cosh —A 


- p 2 ^ 2 l(mo 2 " mg 2 )^ - mo 2 P 2 


apino 


[tain + Mg(d - a)*| 2 , d(l + mQmgP 2 ) - a 

-t =!— cosh —A — 


2(“o + “2) 


P jtaiQ + mg(d “ a) 


fdrriQ + Dip(a - d)1 2 , d(l “ mQmgP 2 ) - a 

- -is — cosh x — 


2(mo - mg) 


P jduiQ + mgCa - d) j 


and 
3 > = 


1 I ]+a 2 m 1 ^m2\l “ mg 2 p 2 ^ am^l - m^P 2 ) - 2am ± 

< - — — . . — cosh — — — — 


Jl - m^p 2 (n^ 2 - “l 2 )^ 1 • m i 2 P 2 


2a pm^ 


[dmo(mi_ - mg ) + 2am 1 mgj 2 (l ~ m^rngP 2 ) - Sam^ 


2mi(mi - mg) 


cosh 


P |dm 0 (m 1 - mg) + 2am 1 mgj 


[tang (mi + mg) - 2am 1 mg _ j 2 _2_ dmo(l + m 1 mgP 2 ) - 2am 1 
“ - cosh — — — 1 " ■ 


+ “2) 


Pjdmg^i mg) - 2am 1 m 2 j 



Wing Drag far Taper Ratio 0 
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For i- < p < — , use A + B + C where 


2a W m l , _„ v -l 1 - ¥/ 

A s ■ 11 " ' ~~ ' " 1 JUq { eu j cosh - 

^1 - ICq 2 ? 2 ^! - IDq)^, 2 - % 2 )_ 


2m 1 (m 0 - m x ) log (l + \jl - 


a^jmo P 

^(^0 “ m l) lo S 

1 - m 1 2 3 2 (m 0 2 - m ± ? - ) ^ 


1 - 

- 2m-]_ cosh 

P( m 0“ m l) 


(l + Jl - m-^p 2 


a ^Hiffln 2 


Jl - m 2 2 P 2 (2m 1 - ihq) (mo - (3m! - ) 


_ 2 _„-l m l( x “ ^ 2 P 2 ) + m l * “0 

2nin cos - ■ 

^(*0 - r ^) 


+ (£sit_ - D3o)(mo “ %) cos 


-1 1 


- - mo) cos 


_! m-^l - m^jtqQp 2 ) + - 

apm^ino - mi) 


For p = — , use A + B plus 
“2 


„ 2 1 V m l / 

a “l^O + 7 - *o 

3% " “0 ( 2 m 1 - m 0 )(n b - m^ 3 ' 2 - afc 
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Far p » — , use B + C plus 
®0 


2a 




2m^ 


(=b “ -i)r - ^x 2 ? 2 (^i * “teX 1 ^ + m i) 


Section Dreg for Taper Eatio 1.0 


c d w c *\jl ” 

l«^(t/c) 2 


for 


0 < y < 


aB b 


1 - 0nfc 


» A + B 


< < 2aafo 

1 - pn^j * 1 - pa^ 


A + B + C 


1 - p*b 


< y < » 


•where 
A » 


k(2j + ajao) - y(l + P% 2 ) + amo 12y 1 + P 2 ^ 2 

= ■■ cosh - — =£• cosh - 

fflQ Pmo(2y + ano) no 2pn^ 


+ 4a cosh 


-1 y(l ~ P^b 2 ) + amp 2(y + amp) ^ x y (l 4 P 2 ^ 2 ) + 2agQ 

Wq 2Pno(y + a^o) 


- 2a cosh 


-1 y(l - PV) + 2an^ 


2aPfflQ 2 


4(2y-au^) ^ 1 y(l+p 2 m 0 2)-ara 0 , ^ y(l- P 2 ^ 2 ) -aa^ 




COHI1 


P®o(2y “ a^o) 


■ - #a cosn 


and 
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C = 2a cosh" 1 . S < 7 ' am °) 003h -l 

8a^ s “o 2Pm 0 (y - an^) / 


Wing Drag for Taper Batio 1.0 


Cp„«4 ~ P 2 ^ 

^)(t/c) 2 


0 < dfflQ < 


A + B 


amo 2amo 

< d%, < 

1 - Pno * 1 - pc 


5= A + B + C 


1 - PIQO 


< duiQ < w 


•where 


- " P 2 ^ 2 ) + a (3 “ P 2 ^ 2 )] -! (l - P 2 ^ 2 ^ 

i - p 2 ^ 2 ° OSh af3ra 0 


-i n 2 2 

1’P Hq 


\|(l - P 2 m 0 2 )d 2 + 2ad-i- a 2 + 1^(21+ a) 2 cosh’ 1 - ■ + - ~ 1 


Pn^(2d + a) 


- Sd^niQ 


_ x 1 + p^ 


an\)ja(3 " P 2 *^ 2 ) + 2d(l - P 2 ^ 2 ) -1 (l - P^m^d + 2a 

■ a — — — - — ■ - * cosh — — — — 

1 - p 2 n^ 2 2aPm 0 


♦ _ g2_ - % (d ♦ a ) 2 oosh' 1 O-^fV^ai 

Jl - p2^2 2Pm 0 (d + a) 
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B 


tmo 


[a ( 3 - PVQ - ta(l - P% 2 )] -1 (1 - ~ a 

1-eV CO Epino 


31110 — J(l - p^nu^d 2 - 2s d + a 2 4 - ibq( 2& - a) 2 

)|l “ P^o 2 


-1 (l+P^®0^)^“ a 

cosh " " 

**>(2d-a) 


and 


- P^) g ) ~ g&Q- ~ P^o g 3 co0h -l (l “ PW) d " 2a 

1 - P% 2 


- , ^ ^j(i - p 2 ^ 2 )* 2 - 4ad + 4a 2 

f * P% 2 


o _T (1 + P 2 ^ 2 ^ - 2a 

- BL(d - a) 2 cosh' 1 ^ 

2PmQ(d - a) 


Vhen p 


i- , use the folloving expression 
®0 



[fed * a)3/ 2 - (d ♦ a)3/ S ] 

3nd\^y. 
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APPENDIX B 


EVALUATION OF EQUATION (7) FOE TIP EFFECTS 


— *f . 

PECTS p < — 

_ -m 0 _ 


Section Drag Increment for 0 < Taper Eatio < 1 


For aspect-ratio limitations, see equations (6) 


C(i tip crt . L djn o 2 ( 1 + P m l) " am i“0 

- — = A I for -i— < y < duu 

Mt/c) 2 \ mi (l + Pkq) ~ 


+ B ( for 


during (1 + Pm 1 ) - 




. „ A™ “M 1 + P»t>) - amg „ „ , 

+ C for < y < dru 

V 1 + Pm 2 - - 


vhere 


A = -2 (y - dn^) cosh 


-1 7% + am^mQ - aniQ £ 
pm 1 m 0 (dm 0 - y) 


+ 2 [yfap ~ m i ) " , 1 ym^l - roQmjp 2 ) - dmo 2 (l - m-^P 2 ) + am^ 


mJl - m/p 2 


Pmi|y(mo - ^ ) - am^j 
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where 


A = 


\|l “ P 2m 0 2 , . , . -i 7 + ®o(2a - d) 

-- -{y - du^j) cosh 

Portae “ y) 


“0 


+ 2a cosh 


-i C 1 ~ ap g P g )(y * dm p) -■ 2am o 

2a pu I q 2 


and 


- P^ 2 • v -i y “ fflo(4 " a) 

B a (y - dn 1 q ) cosh 


“o 


- ha cosh 


P^(dB b - y) 


-1 (l “ mQ 2 P 2 )(y - amp) + amp 
apB^ 2 


There is no tip effect whatsoever for the wing of taper ratio 0, 
and the increment in wing wave drag caused by the tip is identically 
equal to zero for all cases satisfying these aspect -ratio limitations. 
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Case II 


11 2 mo 

— < A < - and A < — 

2p ~ P 1 + Pmo 


c Dtip ^V^ 


^(t/c)' 


Odmo 


Jo 


y-'-amp 

“0 x *+«-a + a 0 P 2 (y +aa to) 

cosh '■ " ■ * — dx ay 

P (y + hjqx + am o ) 

(y+dn|Q)P^d-a 


= 11^(8 - d)‘ 


.]_ a - cta^P 
cosh - cosh 


2 p 2 _ x 2 a - a(l+ mo 2 P 2 ) 


Pidq |a - d | 2 pmo|a - a 


c 

- a jBq cosh 


2a - a(l - m^p 2 ; 

23P131Q 


o _i a + dmo 2 p 2 

+ m_ (a + a) cosh — 

U pm 0 (d + a) 


+ ^2vll - p 2 m 0 2 /cosh" 1 2a --— - 4 cosh" 1 — 

v ^ \ p^a PV. 


2 



and A > 


K'-CA KM Ho. L ; T E23a 


31 




NACA m No. L7E23a 


0 .-.a-dmo^P 2 p -i a - 2dmQ 2 P‘ 

. (a -d) 2 cosh 1 — — •- (a - 2d)^ cosh - — ~ 

T PhIq |a - d | | a - 2d 


, a+dmQ 2 p 2 

+ (d+a)^ cosh 

Pmq( d + a) 


p >1 a + 2dmQ 2 p 2 
- ( 2 d*ar cosh - r 


Pa^(2d + a) 


- 2& 2 \jl - P^tUQ 2 ^cosh 1 * 


a i q -1 & 

— k cosh — 

Pm 0 a 2Pm 0 d 


1 mo 20Q 

0 < A < — and ^ A < 

23 1 + P^o 1 + P 111 © 

The lover limit for y is changed to 0 in the first integral 
of case IH and the resultant expression for the drag is 


W 1 ' p2 "fa ; 

&no(t/c) 2 


o -x a “ amo 2 P 2 


= dvn <(a - d)“ cosh 

1 |_ Pm^a-d 


- cosh 


ml 2a - d(l + mo 2 p 2 ) 
2Pm^ja - dj 


p _■) a ~ 2drn 0 ‘-p 2 « ^a+am^P 2 

- (a- 2d) 2 cosh 1 +(d+a) s cosh - 

PmQ|a-2d Pm^(d+a) 


- (2d 


„ . a + saa^P 2 e -i 2a - d(l - eVO 

+ a) cosh - a cosh 

Pm^(2d+a) Sapm^ 


d 2 \ll - "p^nu 2 ( 2 cosh -1 — 8 cosh 1 

” l Px^d 2pm 0 d 


1 , -1 2a - d 

- cosh 
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Case V 


1 ®0 
0 < A < — and A < — — — — 

2P 1 + PxDq 


All four lover limits for y are changed to 0 in case ill* 


^tin^f ‘ p2 “0 2 J, ,*2 i ,-l a " 

— — — =mQ<(a - d) cosh 

am 0 (t/c) 2 I L pa o! a " d! 


i-, 2a - d(l + mn c 3 2 ) 2 -l a " 2,3ln O P 2 

- 2 cosh 1 — — -- - (a - 2d) 2 cosh A 

2Pmo|a - d I pm^a -2dl 


, a - d(l+ 3 2 idq 2 )^ 5 , a+ dmo 2 P 2 

•cosh x -■ +(d + a.) cosh >' 1 

Pm 0 |a-2d| J Pm^d + a) 


0 a+2dffio 2 P 2 p 

(2d+ a) 2 cosh 1 + a 


Pxa 0 (2d+ a) 


+ a cosh 


-1 3 ‘ d ( X " V^ 2 ) 


apmQ 


2 cosh 


ml 2a - d(l - P% 2 ) 

t 

aapr^ 


- d 2 \jl - p^rru 2 ( 2 cosh 1 - — - 


P®gd 


Pmgd 


■+ 2 cosh 


-1 a 


Pn^d 


2dPmQ 
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Figure 1 .- Symbols and distributions of sinks and sources for a tapered wing. 





value I y-va/ue. 
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Fig. 3 



Figure 3.- Information pertinent to integration limits in equations (5) and (7). 




Figure 4 .- Additional tip effects for wing of constant chord (taper ratio 1 . 0 ) 
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Fig. 5 



Figure 5.- Family of tapered plan forms used for calculations. Plan forms shown have 

same area. 
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Fig. 7 



Figure 7.* Section wave -drag distribution for wing of taper ratio 0. Mach number, 1.414; 

aspect ratio, 8.92; sweepback angle, 60°. 



Fig. 8 
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Figure 8.- Section wave -drag distribution for wing of taper ratio 0.5. 
Mach number, 1.414; aspect ratio, 1.63; sweepback angle, 60°. 
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Fig. 11 
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Figure 11.- Variation of wing wave-drag coefficient with Mach number. Sweepback angle, 50 
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Fig. 15 
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Figure 15.- Variation of wing wave-drag coefficient with taper ratio for three families of wing 

plan forms. Mach number, 1.414; sweepback angle, 60°. 



